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3)  Our  VERB  code  successfully  reproduced  the  behavior  of  the  radiation  belts  during  the 
Halloween  solar  storms  in  2003.  We  have  also  showed  that  during  extremely  strong 
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will  populate  the  inner  belt  (below  1.5  Re)  ,  where  enhanced  relativistic  electron  fluxes 
may  persist  for  several  years  [Shprits  et  al,  2010]. 
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Technical  report 

1.  Bounce  resonance 

We  have  computed  scattering  rates  due  to  ULF  waves,  day-side  and  night-side  chorus 
waves,  EMIC  waves,  hiss  waves  in  the  regions  of  plumes  and  inside  the  plasmasphere, 
lightning-generated  whistlers,  and  whistlers  produced  by  VLF  transmitters. 


a  .  deg  a  ,  deg 

eq  eq 


Figure  1.  Bounce  and  drift  averaged  pitch-angle  diffusion  coefficients  for  scattering 
electrons  at  0.1,  0.5  and  1  MeV,  due  to  resonant  scattering  by  plasmaspheric  hiss, 
lightning-generated  whistlers,  and  VLF  transmitters. 

We  have  also  showed  that  another  scattering  mechanism  (bounce  resonance  scattering) 
driven  by  magnetosonic  waves  plays  a  very  important  role  for  scattering  radiation  belt 
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electrons  [Shprits,  2009].  This  scattering  mechanism  was  originally  proposed  by  Roberts 
and  Schulz  [1968]  for  fast  Alfvc'n  waves  and  has  not  been  studied  since  then. 

Figure  1  (left  row)  shows  that  pitch-angle  scattering  outside  of  the  slot  region,  at 
L= 3.5,  can  provide  rather  weak  but  still  non-zero  scattering  of  electrons  up  to  90°  pitch- 
angle.  In  the  slot  region  at  L= 2.5,  the  combination  of  whistler-mode  waves  cannot  affect 
electrons  above  -88-89°  for  all  considered  energies  of  100  keV,  500  keV  and  1  MeV. 


Figure  2.  Harmonics  of  the  bounce  resonance  frequencies  as  a  function  of  L-shell. 
Colored  bold  lines  show  Hydrogen,  Helium,  and  Oxygen  gyro  frequencies. 

Figure  2  shows  the  first  to  the  twentieth  harmonics  of  the  bounce  frequency  as  a  function 
of  L-shell  computed  for  the  dipole  field.  The  bold  red  line  shows  hydrogen 
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gyrofrequency,  computed  for  the  dipole  magnetic  field.  The  bold  blue  line  shows  the 
lower  hybrid  frequency.  For  a  broad  range  of  energies  from  200  keV  to  5  MeV  and  a 
broad  range  of  L-shells,  multiples  of  the  bounce  frequency  lie  between  the  hydrogen 
gyro-frequency  and  the  lower-hybrid  frequency  where  magnetosonic  waves  are  observed. 
The  absence  of  scattering  rates  for  equatorially  mirroring  electrons  indicates  that  another 
scattering  process  affecting  electrons  mirroring  very  closely  to  the  geomagnetic  equator 
is  required  to  explain  the  global  coherency  of  electron  fluxes.  In  this  study,  we  showed 
that  magnetosonic  waves  and  EMIC  waves  can  violate  the  second  adiabatic  invariant  of 
nearly  equatorially  mirroring  electrons  and  plausibly  provide  pitch-angle  scattering  for  a 
wide  range  of  electron  energies  and  L-shells  inside  and  outside  the  plasmasphere. 
Violation  of  the  second  adiabatic  invariant  will  also  result  in  the  random  change  of 
parallel  energy,  which  may  also  provide  an  in-situ  acceleration  of  electrons  by  means  of 
local  energy  diffusion. 

2.  VERB  code  simulations  of  the  Halloween  storm 

The  computed  diffusion  coefficients  were  incorporated  into  our  3D  Versatile  Electron 
Radiation  Belts  (VERB)  diffusion  code  [Subbotin  and  Shprits,  2009;  Shprits  et  al,.  2009]. 

VERB  diffusion  code  can  reproduce  and  explain  the  dynamical  evolution  of  the 
electron  radiation  belts  during  one  of  the  most  unusual  episodes  of  geomagnetic  activity 
ever  observed  in  space  (Figure  3).  Since  in  its  current  formulation  the  VERB  code  does 
not  use  measurements  for  any  of  the  boundary  conditions  or  initial  conditions  and  is 
driven  only  by  the  geomagnetic  index  Kp,  the  VERB  code  can  predict  the  possible 
consequences  of  an  even  stronger  geomagnetic  storm  that  may  occur  in  the  future  or 
veoccurred  in  the  past. 
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SAMPEX  obsenarions 


10/27  11/01  11/06  11/11  11/16  11/21 


\TRB  simulations 


300  305  310  315  320  325 


Day  of  the  Year,  2003 

Figure  3.  (a)  SAMPEX  observations  of  2  -  6  MeV  electron  fluxes  in  (logio(cm"“  sr'  sec" 
!))  during  the  2003  Halloween  storms,  DOY  300  -  326  (October  27  -  November  22)  1/2 
day  averaged.  Fluxes  are  presented  as  a  function  of  time  and  distance  from  the  Earth  in 
Re  (Earth  radii).  Averaging  is  applied  to  eliminate  periodic  variations  associated  with  the 
spacecraft  passing  over  the  regions  with  weak  magnetic  fields,  usually  referred  to  as  the 
South  Atlantic  Anomaly,  (b)  VERB  code  simulations  of  the  dynamics  of  the  radiation 
belts’  2-  6  MeV  electron  omni-directional  fluxes  modeled  at  the  altitude  of  850  km  in  a 
dipole  field  in  (logio(cm'2  sr"1  sec"1)).  A  one-day  moving  average  is  applied.  Unlike  in  the 
radial  diffusion  studies  that  used  data  at  the  outer  radial  boundary,  the  VERB  code  used 
for  this  study  does  not  utilize  observations  and  is  driven  by  the  Kp  index  only.  The  outer 
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boundary  condition  is  set  up  at  L  =  7  and  its  variation  is  parameterized  by  the  Kp  index, 
(c)  Evolution  of  the  Kp  index  used  for  parameterizations  of  waves  and  modeled 
plasmapause  location  Lpp  (24). 

Our  observations  of  the  radiation  belts  in  space  are  limited  to  only  3-4  solar 
cycles.  Coronal  mass  ejections  (CME)  on  the  Sun  may  produce  even  stronger  storms  than 
those  observed  in  space  over  the  last  3-4  decades.  While  the  in-situ  observational  data  set 
is  a  very  short  indirect  measurement  of  the  intensity  of  a  storm,  taken  from  ground 
measurements  or  visual  measurements  date  back  to  the  17th  century.  While  ground 
disturbance  during  a  regular  storm  is  on  the  scale  of  -50  to  -100  nT,  during  the  recent 
Halloween  super  storms,  when  the  slot  regain  between  the  two  belts  was  populated,  the 
disturbances  on  the  ground  were  on  the  scale  of  -400nT. 

Possibly  the  strongest  ever  observed  geomagnetic  storm  occurred  on  Spetermber, 
1859.  It  was  also  the  first  ever  visually  observed  storms.  This  “perfect  storm”  is  usually 
called  the  ‘Carrington’  storm,  after  one  of  the  first  witnesses  to  observe  and  document 
this  storm  [Carrington,  1895].  Using  measurements  of  the  magnetic  field  component 
from  the  Clabea  Observatory  in  Bombay/Mumbai  taken  during  1846-1867,  Tsurutani  et 
al.  [2003]  inferred  that  disturbances  on  the  ground  produced  by  such  a  storm  may  have 
exceeded  1600  nT,  which  is  approximately  four  times  as  large  as  was  observed  during  the 
recent  Halloween  storms.  Other  strong  storms,  which  were  likely  weaker  than  the 
‘Carringtorn’  storm,  occurred  in  1890,  1909,  1946,  1956  and  1960,  all  before  extensive 
in-situ  observations  were  available  [Shea,  2006]. 
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Figure  4.  (a)  VERB  code  simulations  of  the  dynamics  of  the  radiation  belts  ’2-6  MeV 
electron  omni-directional  fluxes  at  an  altitude  of  850  km  in  (logio(cm"“  sr'  sec'  ))  for  the 
2003  Halloween  storm,  (b)  same  as  (a)  but  during  a  superstorm  with  Kp  used  for 
parameterizations  increased  by  2  at  all  times  when  the  Kp  index  is  greater  than  7.  (c) 
Evolution  of  2  -  6  MeV  radiation  belt  electron  fluxes  at  an  altitude  of  850  km  for  L  =  1 .5 
and  L  =  2  normalized  to  the  maximum  value  of  the  fluxes  for  a  modeled  perfect  storm, 
(d)  Evolution  of  the  driver  of  the  simulations  derived  from  the  Kp  index  and  modeled 
plasmapause  location  Lpp.  The  original  time  series  of  Kp  index  used  for  the  simulations 
described  in  section  2  was  increased  by  two  when  Kp  index  exceeded  7. 
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Figure  4A  shows  the  original  simulations  (same  as  Figure  3B)  presented  here  for 
comparison.  Figure  4B  shows  the  results  of  the  simulations  of  a  hypothetical,  stronger 
than  Halloween  storm  and  compares  them  to  the  original  simulations  (Figure  4A) 
described  in  detail  in  the  previous  section.  While  the  main  phase  drop  out  of  fluxes  in  the 
simulated  ‘perfect  storm’  are  stronger  then  for  the  simulated  Halloween  storm,  electrons 
can  still  penetrate  to  lower  L-shells,  deep  into  the  inner  radiation  zone.  Severe  erosion  of 
the  plasmasphere  would  allow  ULF  waves,  which  may  be  otherwise  reflected  at  the 
plasmapause,  to  penetrate  to  such  low  L-shells.  These  electrons  are  further  accelerated  by 
the  local  acceleration  and  result  in  stronger  2-6  MeV  fluxes  and  a  larger  special  extent  of 
fluxes.  The  modeled  storm  on  November  20  also  diffused  electrons  to  such  low  L-values 
and  further  increased  fluxes  below  L=1.5 

The  main  difference  between  the  evolutions  of  fluxes  is  that  while  at  the  outer  edge 
of  the  inner  zone  fluxes  decay  on  the  time-scale  of  5-10  days,  fluxes  in  the  heart  of  the 
inner  zone  remain  practically  unchanged  during  the  simulation  of  the  October  29th  storm. 
A  modeled  consequence  of  a  new  flare  on  day  324  further  increases  fluxes  at  L=1.5. 
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Figure  5.  Illustration  of  the  evolution  in  the  Earth’s  radiation  belt  fluxes  due  to  a 
superstorm  that  compresses  the  plasmapause  to  1.5  Earth  radii.  Color  indicates  the 
intensity  of  the  radiation  belt  fluxes,  (a)  Typical  two-zone  structure  of  the  Earth’s 
radiation  belts  (b)  acceleration  of  particles  during  the  superstorm,  intensification  of  the 
inner  belt  and  depletion  of  the  outer  belt,  and  (c)  increased  inner  belt  fluxes  after  the 
superstorm  that  may  persist  for  several  years.  The  green  color  indicates  the  average 
amplitude  of  the  radiation  belt  fluxes.  The  blue  line  shows  the  plasmapause,  which  is 
usually  outside  the  slot  region  and  assures  that  local  acceleration  may  only  increase 
fluxes  in  the  outer  radiation  belt.  During  the  Halloween  storms  the  plasmapause  was 
compressed  to  only  2  Earth  radii,  which  did  not  allow  electrons  to  penetrate  into  the  safe 
zone  below  1.5  Earth  radh. 

Fig.  5  illustrates  the  general  structure  of  the  radiation  belts  before  (a),  during  (b), 
and  after  (c)  such  storms.  During  a  superstorm  as  described  above,  the  plasmapause  is 
eroded  and  compressed  below  1.5  Earth  radii  which  allows  ULF  waves  to  penetrate  into 
the  inner  belt  and  diffuse  electrons  radially  much  closer  to  the  Earth.  Local  acceleration 
in  the  inner  belt,  where  the  magnetic  field  is  strong  and  the  plasma  density  is  low, 
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enhances  inner  radiation  belt  fluxes  and  will  make  them  several  orders  of  magnitude 
larger  than  the  outer  radiation  belt  fluxes,  depicted  as  a  red  belt  in  Fig.  3  (b)  and  (c). 

These  results  are  consistent  with  early  observations  of  the  electron  fluxes  after  the 
Starfish  nuclear  detonation,  which  populated  an  artificial  radiation  belt  at  L-1.3  as 
observed  by  Injun  1  [O’Brien  et  ah,  1962].  Observations  of  Paulikas  et  al.  [1967]  yield 
lifetimes  of  ~  265  days  at  L=1.3  for  electrons  above  4.5  MeV.  Such  a  super  storm  may 
significantly  decrease  the  lifetime  of  the  polar-orbiting  and  low-orbiting  satellites  which 
traverse  the  Earth’s  inner  radiation  belts. 


3  Simulations  of  the  Slot  Region  Dynamics 

In  this  section,  we  present  simulations  of  the  dynamics  of  the  slot  region  and  compare 
them  with  observations. 


Steady-state  solution  (Kp=2) 


Figure  6.  Initial  radial  profiles  of  fluxes  for  various  energies  obtained  from  a  solution  of  the 
steady  state  radial  diffusion  equation. 


For  the  initial  condition,  we  used  a  steady- state  solution  of  the  radial  diffusion 
with  Kp=  1  by  considering  both  [I  think  this  is  being  redundant] components  of  the  radial 
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diffusion  coefficients  with  the  simplified  lifetimes  described  in  section  2.1,  which  is 
displayed  in  Figure  6  for  various  energies  as  a  function  of  L*.  In  Figure  6,  one  can  see  the 
tendency  that  for  a  higher  energy  the  slot  region  is  more  pronounced  and  is  located  more 
inwardly  as  already  mentioned.  Based  on  the  initial  flux  profiles,  we  have  also  solved 
both  the  1-D  and  2-D  diffusion  equations  and  compared  them  to  each  other. 

Based  on  the  assumptions  above  and  by  solving  the  Fokker-Planck  equation  for 
the  radiation  belts,  we  have  conducted  a  2D  diffusion  simulation  of  the  time  evolution  of 
the  IMeV  electron  flux  for  an  equatorial  pitch  angle  ao~90°.  The  result  is  displayed  in 
Figure  7c.  For  comparison,  we  also  present  the  corresponding  CRRES  observations 
(Figure  7a)  and  a  1-D  radial  diffusion  simulation  (Figure  7b).  In  Figure  7c,  one  can  see 
that  our  diffusion  model  predicts  the  instantaneous  location  of  the  upper  boundary  of  the 
slot  region,  a  slot  region,  and  stable  inner  belts,  which  shows  a  general  good  agreement 
with  the  CRRES  observations. 

Now  we  model  more  realistic  situations  by  taking  into  account  several  factors: 

(1)  the  motion  of  the  plasmapause  in  response  to  the  Kp  variation  (2)  the  flux  variation  at 
the  outer  boundary  (L*= 5.5)  as  shown  in  Figure  5  (middle)  and  (3)  pitch-angle  scattering 
by  day-night  chorus  waves  outside  the  plasmapause  for  which  the  bounce-averaged  pitch 
angle  diffusion  coefficients  are  calculated.  The  readers  are  referred  to  the  work  by  Shprits 
et  al.  [2009]  for  more  detailed  descriptions  of  the  chorus  model.  For  comparison,  we  have 
also  conducted  a  1-D  radial  diffusion  simulation  with  the  empirical  lifetimes  taken  from 
Shprits  et  al.  [2006] . 
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Figure  7.  Time  evolution  of  IMeV  electron  fluxes  with  equatorial  pitch  angle  of  «0~90°  based 
on  (a)  CRRES  MEV  IMeV  observations  (b)  a  1-D  and  (c)a  2-D  diffusion  simulation  with  a 
variable  outer  boundary. 
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